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Abstract. In this paper we derive stellar parameters 
for the Wolf-Rayet star in the 7 Velorum binary system 
(WR11), from a detailed non-LTE model of its optical and 
infrared spectra. Compared to the study of Schaerer et al., 
the parameters of the WC8 star are revised to a hotter ef- 
fective temperature (T e g ~ 57 kK), a higher luminosity 
(log(L/L Q ) = 5.00), and a lower mass loss rate (log(M 
/ M Q /yr) = -5.0 - using a 10% clumping filling factor. 
These changes lead to a significant decrease in wind effi- 
ciency number, from 144 to 7, so that the driving mecha- 
nism of the wind of this WR star may be simply radiation 
pressure on lines. The derived spectroscopic luminosity is 
found to be 40% lower than that derived by De Marco 
& Schmutz through the mass-luminosity relationship for 
WR stars (log(L/L ) = 5.2). 

The paper furthermore presents a comparison of the 
independently-developed modelling programs, CMFGEN 
and ISA- wind. Overall, there seems to be very reasonable 
agreement between the derived parameters for WR11, ex- 
cept for the carbon content, which is 2 times higher for CM- 
FGEN (C/He=0.15 vs. 0.06, by number). The comparison 
also confirms a disparity in the predicted flux at A<400 
A, found by Crowther et al., which will have effects on 
several nebular line strengths. 

The paper also presents the first independent check of 
the photon loss mechanism proposed by Schmutz. We con- 
clude that, not only is it important to include very many 
lines to realistically model line blanketing, but in particu- 
lar those ones that critically interact with strong resonance 
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lines (e.g. He 11 A303.78). The inclusion of these latter lines 
may significantly alter the wind ionization structure. 
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1. Introduction 

The double-lined spectroscopic binary 7 Velorum (WC8 
+ 07.5) can be used to derive the mass of the Wolf-Rayet 
star, a quantity which cannot be directly determined in 
the case of single objects, but which provides an important 
constraint in computations of the late stages of massive 
star evolution. From it, and the mass- luminosity relation- 
ship for WR stars (Schaerer & Maeder 1992) a value for 
the WR luminosity can be derived and compared to that 
derived from spectroscopic analyses. 

Orbital parameters together with the mass ratio were 
determined by Schmutz et al. (1997), while the 7.3a HIP- 
PARCOS distance allowed Schaerer et al. (1997) to derive 
the system's total mass (.M(0+WR)=29.5±15.9 M Q ), the 
absolute visual magnitude (My (0+WR)=— 5.39) and to 
carry out a preliminary spectroscopic analysis of both 
components. Schaerer et al. (1997), however, modelled 
only the WR Hei A4471 / Hen A4541 line strength ra- 
tio, assuming the mass-loss rate and the carbon abun- 
dance. Additionally, line blanketing was not accounted for 
by their model and their calculation used simple helium 
and carbon model atoms. Their stellar parameter deter- 
mination should therefore be considered preliminary. 

De Marco & Schmutz (1999; hereafter Paper I) car- 
ried out a full spectroscopic analysis of the O star through 
which its parameters were determined simultaneously with 
the light ratio between the WR and O stellar components. 
From a solution of the hydrodynamic wind equations, the 



2 



O. De Marco et. al.: The 7 Velorum binary system. II 



mass of the O star was also derived and found to be con- 
sistent with the evolutionary O star mass. From it and 
the mass ratio of Schmutz et al. (1997), the mass of the 
WR star could be calculated (M=9.0±1.0 M Q ). The com- 
bined mass of the system was therefore derived to be 
39 M©, higher than the value of Schaerer et al. (1997), 
but consistent within their uncertainty. From the mass of 
the WR star and the theoretical mass-luminosity relation 
of Schaerer &; Maeder (1992), a value for the WR bolo- 
metric luminosity was derived (log(L/L©) = 5.2), signif- 
icantly higher than the spectroscopic luminosity derived 
by Schaerer et al. (1997; log(L/L©) = 4.8). Finally, using 
the synthetic O star spectrum, the WR11 spectrum was 
de-convolved from the O star component allowing the WR 
emission lines to be measured with greater accuracy. 

In this paper we determine the stellar parameters of 
the WR star using the clumped, line-blanketed stellar at- 
mosphere model CMFGEN (Hillier 1987, Hillier & Miller 
1998). Ultimately, higher spectroscopic luminosities would 
lower the wind performance number, indicating that a 
smaller number of photon scatterings is needed to achieve 
the radiative driving of the wind. This would facilitate the 
derivation of the wind velocity law. 7 Velorum provides us 
with a new tool to test the model luminosity with an in- 
dependent method in such a way as to impose a further 
check on the wind momentun number. 

A second aim of this paper is to continue the compari- 
son of the non-LTE stellar atmosphere code CMFGEN and 
the independently developed Sobolev approximation code 
ISA- WIND (de Koter et al. 1993, 1997), which is used in 
combination with the line-blanketing Monte Carlo code of 
Schmutz (1991). This study follows from the comparison 
carried out by Crowther et al. (1999) for the WN8 star 
WR124. 

Finally, Schmutz (1997) showed that line-line interac- 
tion between He 11 Lya and metal lines at similar wave- 
lengths can lead to a larger spectral luminosity. The in- 
clusion of this interaction in the models, an effect that 
he called "photon loss", showed that WR6 might have 
a larger luminosity than previously thought, and allowed 
him to calculate a velocity law for its outer wind. Since 
line-line interaction is included in the CMFGEN model code, 
we can test whether the approximation of Schmutz (1997) 
is valid and quantify the importance of the photon loss 
effect in the case of WR11. 

In Sect. we summarise the observations, while in 
Sect. I we discuss the light ratio between the O and WR 
stars and the reddening towards the system. In Sect. |J we 
describe the model atmosphere codes used in the determi- 
nation of the stellar parameters. We present our results in 
Sect. together with the code comparison, while in Sect. || 
the effect of photon loss is studied. We finally draw our 
conclusions in Sect. 



2. Observations and reduction 

Our optical observations of 7 Velorum are the same as 
those used by Schmutz et al. (1997), where a full account 
of the data reduction can be found. In summary, several 
spectra have been obtained at the 50 cm ESO telescope 
in conjunction with the Heidelberg Extended Range Opti- 
cal Spectrograph (heros) in the ranges 3500-5500 A and 
5800-8600 A covering the binary period. The resolution 
is R=20000, while the S/N ratio ranges between 100 at 
3600 A and 250 at 6000 A. The phase-averaged spectrum 
of the 7 Velorum system was rectified as described in Pa- 
per I. It is important to remind that rectification of WC 
stars is a very delicate operation because of blending of 
large numbers of broad lines, which leave few identifiable 
continuum points. Moreover, spectral ranges at the end 
of instrumental orders (5300-5500 A and 5800-5900 A), 
suffer from additional uncertainty. Unfortunately this af- 
fects the diagnostic lines Hei A5876, Hen A5412, and Civ 
A5471. Extreme care was taken in rectifying these ranges 
and consistency between these lines and those residing at 
other wavelengths was insured. 

Observations of the 10830 A Hei line were obtained at 
the European Southern Observatory (ESO) New Technol- 
ogy Telescope on January 12 and 13, 1996, with the emmi 
instrument, emmi was used in the remd mode with grat- 
ing No. 7, a slit-width of 1.5 and the spectrum recorded 
on ESO CCD No. 36. Two exposures, each lasting 10 s, 
were taken. Wavelength calibration was obtained with re- 
spect to a ThAr arc lamp, to achieve a resolution of 2 A. 
Relative flux calibration was obtained with respect to a 
180 s exposure of the standard star 9 Crt. 

Mid-IR data of WR11 were obtained as part of 
Guest Observer programme WRSTARS (P.I. van der 
Hucht), with the Short Wavelength Spectrograph (SWS; 
de Graauw et al. 1996) on board the ESA Infrared Space 
Observatory (ISO; Kessler et al. 1996). The SWS AOT6 
observing mode was used to achieve full grating resolu- 
tion, A/AA ~ 1300 — 2500, with the wavelength coverage 
2.38-45.0/im. The observations were flux calibrated to an 
accuracy of 5% and ~20% for the low and high wave- 
length ends, respectively. The possibility that the nearby 
7 1 Velorum might have contaminated the ISO observa- 
tions of 7 Velorum was investigated, but with an aperture 
of 20" centered on 7 Velorum it is not possible that light 
from 7 1 Velorum, at 40" distance, could have entered the 
aperture. For full details on the observations and data re- 
duction procedures, we refer the reader to Morris et al. 
(2000). 

Finally, we have obtained a flux-calibrated UV dataset 
for WR11 as follows. WR11 was too bright for absolutely 
flux calibrated observations with IUE/LORES. Conse- 
quently we have matched IUE observations from July 
1988, to low resolution S2-68 observations from 1973 at 
phase 0.3 (Willis & Wilson 1976). 
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Fig. 1. The light ratio, L\(0)/L\(WR), as a function of 
wavelength. The thinner lines show the error bars derived 
from the original error limits of Paper I 

3. Light ratio correction and reddening 

In this section we discuss the O-WR light ratio of 7 Velo- 
rum. We also re-evaluate the interstellar reddening that 
was adopted in Paper I, namely E(b — v) = 0.03 mag, ob- 
tained from assumed intrinsic colours (van der Hucht et 
al. 1988). 

The spectrum of 7 Velorum used here is the same 
phase-averaged rectified spectrum used in Paper I. Cor- 
rection for the continuum as well as the line contribution 
of the O star was carried out in that paper using the model 

star spectrum and the deduced value of the light ratio. 
A light ratio of L(O) / L(0+WR) = 0.795 ±0.020 was ob- 
tained in Paper I from an average of three lines, namely 
He 1 A4471, Hei A4922 and He 11 A4541. Assuming that 
this light ratio is representative for this spectral range, 
we derive the wavelength dependence of the light ratio 
using synthetic spectra for the O star from Paper I, and 
for the WR star from the following section. The ratio of 
these continuum distributions, adjusted such that L\(0) 

1 La(WR) = 3.8S±g;| at A= 4700 A, is plotted in Fig. |, 
together with error bars. A calibration curve La (O+WR) 
/ La(WR) is created and applied to the optical spectrum 
of WR11 to account for the contribution of the O star con- 
tinuum to the WR equivalent widths. In Fig. || we show 
the effect that the error bars from the light ratio have on 
Hen A5412 and Civ A5471. 

After the spectrum is rectified and light-ratio cor- 
rected, a model can be calculated to fit the WR lines, 
as discussed in the following section. First, we preempt 
the result of this modelling in order to verify our choice of 
reddening. Published broad-band Johnson V magnitudes 
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Fig. 2. Effect of the light ratio uncertainty on the strength 
of important diagnostic lines at He 11 A5412 and C iv 
A5471. (This does not include errors on the synthetic at- 
mospheres used in the construction of the wavelength- 
dependent ratio). The dashed and dotted curves corre- 
spond to the upper and lower error bars, respectively 

for 7 Velorum span a wide range, from V = 1.70 (Cousins 
1972) to V = 1.83 (Johnson et al. 1966). We suspect that 
the broad band observations suffer from contaminations 
by nearby spectral lines such as C iv A5806 and that the 
measured brightness depends crucially on the exact sensi- 
tivity curve of the used setup. Therefore, in order to min- 
imise the contamination of WR line emission, we prefer 
to adopt the narrow-band ubv Smith (1968) photometry. 
Following the corrections by Schmutz & Vacca (1991) to 
the Smith datasets, we obtain v — 1.70 mag and b — v = 
—0.32 mag. Our theoretical combined O+WR model re- 
veals (b — v)o — —0.33 mag, implying E(b — v) = 0.01 or 
E(B - V) = 0.01 mag. 

We have also carried out an alternative reddening de- 
termination, by fitting the combined WR+O model to 
UV (IUE/LORES) and IR (ISO) spectrophotometry, de- 
reddened following Seaton (1979). Fig. || compares ob- 
served and de-reddened spectrophotometry (solid lines) 
with theoretical models (dotted lines), revealing E{b—v) = 
0.06 or E(B — V) — 0.07 mag. Unfortunately, we are un- 
able to reconcile these discrepant reddening values, and 
so adopt E(b - v) = 0.035 mag (E(B - V) = 0.04), in 
good agreement with van der Hucht et al. (1988). Note 
that for this reddening and the adopted light ratio, the 
combined O+WR model flux lies ~10% higher than the 
UV IUE observed fluxes and a similar value below the IR 
ISO observations. 

The distance modulus of 7.06 mag (corresponding to 
a distance of [258^^] pc, hipparcosQ and de-reddened 

The HIPPARCOS distance to 7Velorum has been recently put 
in doubt (Pozzo et al. 2000) by the finding of a new population 
of pre-main sequence stars which might be associated with the 
binary. This new value of the distance re-positions 7 Velorum 
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Table 1. O and WR and system parameters, updated 
from Paper I. From these the mass and luminosity of the 
WR star are determined. Errors on the distance are not 
taken into account - see text 



AMa4700 


1.47± 0.13 mag 


/A4700(O)//A4700(O+WR) 


0.795 ± 0.020 


Mv(0+WR) 


-5.5 ± 0.3 mag 


i 


63±3 deg 


M v (0) 


-5.10 ±0.1 mag 


Teff(O) 


35 000±300 K 


N[He](0) 


0.087 (by number) 


R(0) 


12.4±1.7 Rq 


L(0) 


(2.1±0.3)xl0 5 L 


M(0) 


30±2 M 


M(O) 


(1.78±0.37)xl0~ 7 Mgyr" 1 


^oo(O) 


2500±250 km s" 1 


age(O) 


(3.59±0.16)xl0 6 yr 


My(WR) 


-3.76 ± 0.2 mag 


M(WR) 


9.5±1.0 M Q 


L(WR) 


(1.7±0.4)xl0 5 L 



-a) 



':>t 

0) 



-9 



-10 



-11 



O 



^ -12 



43 











0:WR=3 


8:1 at 4700A - 




E(b-v) = 


= 0.06 


- IUE/L0RES 













3.0 3.5 4.0 4.5 

Log Wavelength (A) 



5.0 



v = 1.56 mag imply A/„(0+WR) = -5.50. The O/WR 
continuum light ratio is 3.61 at the w-band (Fig. [y) so 
that M„(0)=-5.23 and M„(WR) = -3.84. We may con- 
vert the Smith (1968) v to Johnson V via V = w-0.365(fe- 
v) + 0.007, adapted from Turner (1982) to take account of 
the ubv re-calibration by Schmutz & Vacca (1991). From 
theoretical continuum O and WR distributions we obtain 
(b — v) = —0.34 and (b — v) = —0.20, respectively, re- 
vealing absolute line free V-magnitudes listed in Table |l|. 
Relative to Paper I, the difference in the absolute magni- 
tude of the O star, My(O), is —0.04 magnitude. The errors 
listed in Table [l] are solely from the reddening and light 
ratio uncertainties. Errors in the distance, would add an 
additional 0.3 mag in the My uncertainties. Propagating 
these errors onto the derived magnitudes and luminosities 
in Table [l], and subsequently onto the spectroscopic prop- 
erties in Sect. 0, introduces an additional systematic error 
of log (L/Lq)—±0.1 for both cases. In this way we lose 
sight of the actual accuracy reached by this study. 

To be consistent with the slight change in the O and 
WR star V magnitudes we have revised the uncertainties 
on i, and i, while the values of the luminosity and 
mass of the WR star are also slightly modified. 



within the Vela OB2 association at a distance of 360-490 pc, 
in agreement with earlier estimates. Clearly, adopting a larger 
distance would have a dramatic effect on distance dependent 
parameters, such as luminosity, radius and mass-loss (for a dis- 
tance of 360-490 pc, Alog(L/L o )=0.3-0.6). We have, however, 
maintained the 7.3cr HIPPARCOS distance measurement, which 
we trust to be the best distance estimate to date. 



Fig. 3. Comparison of the combined O+WR synthetic 
model (dotted lines) and the IUE/LORES and ISO/SWS 
observations (solid lines) for reddenings of E(b — v) = 0.00 
and 0.06 mag, with a light ratio La47oo(0)/La47oo(WR) 
= 3.8 

4. Stellar atmosphere codes 

In this section we present a summary of the description 
of both codes used (a full account can be found in Hillier, 
1989 and Hillier & Miller, 1998 for CMFGEN and in de 
Koter et al., 1993, 1997, for isa-wind). Along with the 
description of the basic codes, we give an account of the 
Monte Carlo line blanketing code used in conjunction with 
the un-blanketed isa-wind to calculate the additional 
opacity due to the effect of bound-bound transitions on 
the continuum photons. Last, the implementation of the 
photon loss mechanism (Schmutz 1997) into the Sobolev 
code is presented here, for the first time, and a parallel is 
drawn with the potential effect of line-line interaction in 
the line-blanketed code CMFGEN. 

4-1. CMFGEN 

CMFGEN (Hillier 1987, 1990; Hillier & Miller 1998, 1999) 
solves the transfer equation in the co-moving frame, 
subject to statistical and radiative equilibrium, assum- 
ing an expanding, spherically-symmetric, homogeneous or 
clumped, time-independent atmosphere. Line blanketing 
is treated correctly in the transfer problem, except that a 
simplifying 'super level' approach is used (Anderson 1989), 
by which several levels of similar energies and properties 
are treated as a single 'super level'. The population of an 
individual atomic level in the full model atom is deter- 
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mined by assuming that it has the same departure coeffi- 
cient as the corresponding super level to which it belongs. 
For the specific models dealt with in this paper the atomic 
model is shown in Table ||. Elemental abundances other 
than hydrogen, helium, carbon and oxygen, are fixed at 
their solar mass fraction values. 

The stellar radius (R*) is defined as the inner bound- 
ary of the model atmosphere and is located at Rosseland 
optical depth of ~20 with the stellar temperature (T*) 
defined by the usual Stefan-Boltzmann relation. Similarly, 
the effective temperature (T e s) relates to the radius (R2/3) 
at which the Rosseland optical depth equals 2/3. There 
is now overwhelming evidence for the clumped nature of 
WR stars (e.g. Moffat 1999), so we have adopted a sim- 
ple filling factor approach. Following the method proposed 
by Schmutz (1995), we assume that the wind is clumped 
with a volume filling factor, /, and that there is no inter 
clump material. Since radiation instabilities are not ex- 
pected to be important in the inner wind we parametrise 
the filling factor so that it approaches unity at small ve- 
locities. Clumped and non-clumped spectra are very sim- 
ilar, except that line profiles are slightly narrower with 
weaker electron scattering wings in the former. Although 
non-clumped models can be easily rejected, we are unable 
to determine the clumping factor because of the severe 
line blending in WC winds. We can therefore derive only 
M/y/J with our spectroscopic analysis. 

Unfortunately, individual line blanketed co-moving 
frame calculations are generally computationally demand- 
ing, despite the use of super levels, so that a large param- 
eter space cannot be quickly explored. One solution is to 
(i) solve the transfer problem in the classical Sobolev ap- 
proximation rather than the co-moving frame, which the 
code ISA-wind does (de Koter et al. 1993, 1997), and (ii) 
consider line blanketing via Monte Carlo sampling follow- 
ing Schmutz (1991, 1997) allowing the opacity of a huge 
number of lines to be considered. 

4-2. IS A- wind and the Monte Carlo line-blanketing code 

The improved Sobolev approximation code (ISA-wind) 
is described in detail by de Koter et al. (1993, 1997). 
The principal differences with CMFGEN relate to (a) the 
treatment of the line radiation transfer; (b) wind electron 
temperature which assumes a grey atmosphere in local 
thermodynamic equilibrium (LTE); (c) the specific atomic 
model treated, as listed in Table g. 

The co-moving frame method consistently treats a pos- 
sible change of properties of the medium inside the region 
in which line photons can be absorbed and re-emitted. In 
the Sobolev approximation one assumes this line interac- 
tion region is infinitely small. More exactly, the Sobolev 
approximation assumes that the source function and the 
opacities are linear functions such that they can be taken 
out of the transfer integral. This is generally a valid as- 
sumption for Wolf-Rayet stars, where the large velocity 



gradient of the flow results in a small interaction zone. 
The Sobolev classical approximation (Castor 1970) im- 
plies that interactions of line photons with the continuum 
inside of the (in reality extended) interaction region are 
ignored. However, the improved Sobolev approximation 
code ISA-wind does include such continuum interactions. 
Note that use of the Sobolev approximation introduces 
great simplifications in the radiative transfer, resulting in 
an overall iteration process up to about 20 times faster 
compared to the co-moving frame approach. 

Turning to line blanketing, an iterative technique in- 
cluding the Monte Carlo method of Schmutz (1991) is em- 
ployed. The method allows the computation of intensity- 
weighed effective opacity factors, which account for the 
presence of tens of thousands of spectral lines, dominated 
by Fe and Ni. Based on ISA-WIND atmosphere calculations, 
the Monte Carlo code determines the line blanketing fac- 
tors. An iterative procedure is used, such that blanketing 
factors are used by the non-LTE code to calculate a new 
atmosphere, which in turn is used to calculate new blan- 
keting factors. A few iterations are generally sufficient. 
This is due to the fact that the scattering and absorption 
factors are not very sensitive to the exact model intensity. 

The Monte Carlo method deals with the photon scat- 
tering in the correct manner, although no branching is 
taken into consideration, an effect that would tend to re- 
distribute UV photons to longer wavelengths. The ioniza- 
tion equilibrium of metal species is also approximate, in 
that it is determined in relation to the non-LTE structure 
of helium, carbon and oxygen and not by the full solution 
of the statistical equilibrium equations for that species. 

4-3. Photon loss 

"Photon loss" is the name given by Schmutz (1997) to 
the interaction between the photon field of the He II Lya 
line at 303.78 A and nearby metal lines (e.g. the Oiii 
Bowen lines at 305.72 and 303.65 A). This was suggested 
by Schmutz (1997) to be responsible for a decrease in the 
ionization equilibrium of the wind of WR6 (WN4), with 
respect to the equilibrium obtained when no such line-line 
interaction was accounted for. In order to fit the observed 
spectrum when photon loss is included in the model, the 
model's stellar temperature has to be higher, which results 
in a higher spectroscopic luminosity. This contributes to 
a lower wind performance number and allowed Schmutz 
(1997) to calculate a velocity structure for the outer wind 
of WR6. 

Since line- line interaction was not accounted for in any 
atmospheric wind code before Hillier & Miller (1998) im- 
plemented it into the Hillier (1990) code, Schmutz (1997) 
calculated the photon loss factor (i.e. the amount of in- 
teraction between the He 11 Lya and nearby metal lines) 
by comparing the opacity of metal lines in the 303-A re- 
gion to that of the Hell Lya line (his Appendix A). Once 
the amount of interaction is known for a particular model 
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atmosphere, a corresponding fraction of photons can be 
taken out of the radiation field at that particular wave- 
length for all grid points. In this way the same fraction of 
photons is removed from the radiation field at all depth 
points, although in reality the photon loss factor depends 
on the relative opacity of the He 11 Lya line and that of the 
metal lines, which in turn depends on the depth-dependent 
populations. This problem however is not critical as long 
as the opacity of metal lines behaves similarly to that of 
helium. Additionally, the loss of photons is only critical for 
the point in the atmosphere where the wind is recombin- 
ing. We implemented the photon loss approximation into 
the ISA-wind code of de Koter et al. (1993 - see Sect. 4.3.1) 
and used it in the calculation of a model for WR11. 

CMFGEN includes line-line interaction and as such it 
should naturally include the photon loss effect, provided 
the specific metal lines in the He 11 Lya region are included 
in the model atom. We can therefore test whether this 
effect and its consequences for the stellar luminosity has 
been accounted for properly by our approximation. 

4.3.1. Implementation of Photon loss in ISA-wind 

Once the photon loss factor, /, is calculated for a cer- 
tain model atmosphere (see Schmutz 1997, Appendix A), 
a fraction / of the photons is removed from the radiation 
field in the wavelength range 300-304A. To do so, one 
has to modify the expression where the mean intensity for 
this line wavelength is calculated at every code iteration. 
In ISA- wind this was achieved via a modification of the 
Sobolev escape probability. 

In the Sobolev approximation the radiation field J at 
radius point r is related to the line source function S in 
the following way: 



J(r) = (1 - P)S(r) 



(1) 



where (3 is the line escape probability by direct flight and 
where 



S(r) = 



N u A ul 



NiB u i - N U B. 



a I 



(2) 



with N u and Ni the upper and lower level populations, 
respectively, A u i and B u i the Einstein coefficients. Due to 
photon loss, the field at frequencies close to the He 11 Lya 
is modified in the following way: 



J(r)-(l-/)J(r). 

Substituting (1) into (3) we obtain: 

J(r) = (l-(3-f + ff3)S(r) 



(3) 



(4) 



Comparing (5) to (1) we see that the classical Sobolev 
escape probability is modified in the following way: 

/3^(/3 + /-//?). 



5. Stellar analysis 

In this section we present the stellar parameters for WR11 
derived with the two codes CMFGEN and ISA-wind (sup- 
plemented by the Monte Carlo line-blanketing code). We 
will discuss the quality of the fits and the overall flux dis- 
tribution of the model atmospheres. 



5.1. Analysis technique 

In our approach, diagnostic optical lines of Hei (AA5876, 
10830) and Hell (AA4686,5412) are chosen to derive the 
stellar temperature and mass-loss rate. Initially we assume 
the carbon and oxygen abundances to be the same as those 
derived for the single WC8 star WR135 by Dessart et al. 
(2000), because of the similarity of the two spectra (a com- 
parison of the rectified spectra is shown in FigJ|). Once a 
fair agreement is reached between the modelled and the 
observed helium lines, carbon and oxygen abundances can 
be altered to improve the fits to their diagnostic lines. For 
carbon, lines of Cm at 6741 A and of Civ at 5471 A 
are used. For oxygen, the lack of good diagnostic lines in 
our spectral range and the weak and blended nature of 
those that can be identified, make the determination of 
the oxygen abundance impossible. We therefore adopt the 
theoretical O/C number ratio of 0.2 (Meynet et al. 1994). 
Changing oxygen and carbon abundances acts on the to- 
tal opacity of the wind and therefore feeds back on the 
modelled strengths of the helium lines. Parameters have 
therefore to be adjusted so as to retrieve the fits of the 
helium lines. 

We adopt the wind terminal velocity determined by 
St. Louis et al. (1993), — 1550 km s _1 , from the vari- 
ability of the UV spectrum. The velocity law used by the 
CMFGEN code (Hillier & Miller 1999) is a two component 
"beta" law which produces a slower flow than the usual j3 
= 1 velocity law, at points in the atmosphere intermediate 
between the outer regions where v — ► Voo and the photo- 
sphere. Hillier & Miller (1999) found that it leads to better 
fits to the line profiles, while at the same time Schmutz 
(1997) determined a similar velocity law from his hydro- 
dynamic calculation (which, however was carried out only 
for the outer layers of the wind). The two component ve- 
locity law used for our calculations uses vq — 100 km s _1 
(the photospheric velocity), w C orc = 1-0 km s" 1 (the veloc- 
ity at i?*), w xt = 1100 km s _1 (the intermediate velocity) 
and Woo = 1550 km s _1 (the terminal velocity). The two 
beta exponents used are 1.0 and 50, for j3\ and fa, respec- 
tively. The ISA- WIND calculation was performed with the 
same velocity law as CMFGEN. A clumping filling factor of 
0.1 was adopted for CMFGEN. Usually this is derived to- 
gether with the other parameters, by fitting line wings; for 
WC stars the heavy blending makes this task harder. We 
therefore adopted the value obtained by Hillier & Miller 
(1999) and confirmed by Dessart et al. (2000). 
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The atomic data models used for the two codes are 
summarised in Table [2] (see Dessart et al. (2000) for a de- 
scription of the atomic data used in CMFGEN, while ISA- 
wind uses the latest Opacity Project data). Note that 
the number of individual and super levels (Np and Ns) 
mean different things for the two codes: CMFGEN bundles 
Np levels into Ns super-levels, while all ISA-WIND atomic 
levels, listed under the Np heading, are treated individu- 
ally (apart for Hei and Civ for which the upper 32 and 
6 levels, respectively, are super-levels). The data used by 
CMFGEN is more extensive than those used by ISA-WIND. 
Generally this can have an effect on the ionization equilib- 
rium (and hence on the spectrum), or only have an effect 
on particular lines between high lying levels. Tests using 
CMFGEN with reduced atomic data show that while indi- 
vidual lines can be very sensitive to the amount of levels 
used (e.g. C in A8500) the overall ionization of the wind is 
not critically dependent on the size of the atomic model 
provided a minimum amount of levels is used. 

5. 2. Optical and Infra-red line fits 

The overall fit quality of the helium diagnostic lines is fair 
and comparable for the two codes (Fig. ||) . Due to the dif- 
ficulties in rectifying WC spectra and the weakness of the 
lines in the region 4400-5450 A region, it was decided that 
the He I AA5876, 10830 would be better diagnostic lines 
than He 1 A4471, although the fit to this line is not incon- 
sistent. A5876 is very well reproduced by CMFGEN, while 
A10830 is under-predicted by about 15% (one should take 
into account that this line was not corrected for the O 
star absorption line, which can be observed as a depres- 
sion on its blue side). ISA- wind only reproduces A5876 
(albeit not perfectly), due to the sensitivity of A10830 to 
the outer wind electron temperature, which is calculated 
in the grey LTE approximation and fixed for outer radii 
to 10 000 K. 

Hen A5412 is well reproduced by both codes, although 
CMFGEN slightly over-predicts it. Once again we use the 
4440-4600 A region only as a secondary diagnostic: CM- 
FGEN does reproduce the Hen A4541-Civ A4553 rather 
well, while ISA- WIND does not predict the C IV component 
in the blend with an overall poorer result. A4686 is slightly 
under- predicted by CMFGEN, with the ISA- wind synthetic 
line fit being affected by blending with the over-predicted 
C in-C iv line at 4660 A. 

The Hen lines at A4859 and A6560 A are underes- 
timated by both CMFGENand ISA- wind, by about 20% 
and 80% respectively. In the case of A6560 the presence of 
an under-predicted flat-topped C 11 line could be respon- 
sible (see below). It is in fact not uncommon that models 
under-predict line emission from low ionization stages (cf. 
CnA4267), showing that the stratification of the wind is 
imperfectly reproduced. It is possible that this problem 
will be solved once more realistic clumping is included in 
the models. On the other hand there is no blend in A4859 



that could be blamed for the poor fit. For that line there 
could be a slight rectification problem, although this alone 
could not account for the amount by which the model 
undcr-fits the data. 

In Fig. |E| we show a comparison of the observed spectra 
of WR11 and WR135. The two spectra are very similar 
with WR135 having only slightly stronger Hen AA4541, 
4686, 5412 lines and comparable Hei AA5876,7066 lines. 
On the other hand, WR11 has stronger He 11 AA4339, 4859, 
6560, indicating an opposite trend to the rest of the He n 
spectrum! (Hen A4100 is also shown in FigJ| although 
due to the heavy blending with carbon lines, it is difficult 
to make a clear comparison). This points to an anomaly 
in the strength of these Hen Pickering lines. We should 
also note that the ratio of Hen5411 to H/3 is different for 
the two stars, indicating something anomalous in the line 
formation in one or both stars. 

For carbon lines, CMFGEN reproduces strengths and 
shapes to a higher degree of accuracy. ISA-wind synthetic 
lines underestimate C iv AA4440,4780 by about 20% and 
40%, respectively, while Cm AA6741,8500 are weaker than 
the observations by 80% and 50%, respectively. On the 
other hand while CMFGEN reproduces AA4440,4780,6741, 
the Cm line at 8500 A is over-estimated by more than 
100%. The A4660 Cm - Civ blend is overestimated by 
ISA-wind (—60%), while it is underestimated by CMFGEN 
(-10%). Some lines of Cm (e.g. AA8196,8500) and Civ 
(A7063), are under-predicted or not predicted at all by 
ISA-wind partly because they are not accounted for in 
the atomic model. The carbon abundance determination, 
however, rests critically on the ratio of Hen A5412 / C iv 
A5471 which is well matched by both codes. C 11 is under- 
predicted by CMFGEN and not predicted at all by ISA- wind 
(which includes only 7 levels in the atomic model of Cn). 
However the presence of the C II line at 4267 A indicates 
that C 11 lines are present in the spectrum. The spectral 
classification lines C III A5696 and C iv A5806 are not in- 
cluded in the observations, so in Fig. ^ we only show the 
prediction for both models. These two lines are remark- 
ably different in the two models: the spectral types de- 
duced from their equivalent width ratio are WC6 and WC9 
for the ISA-wind and CMFGEN codes, respectively, where 
the criteria of Crowther et al. (1998) have been adopted. 
The model for WR135 of Dessart et al. (2000), with simi- 
lar parameters, has a much larger Civ A5806/Cm A5696 
line ratio, consistent with a WC8 classification. Although 
these two lines form an excellent spectral type diagnos- 
tic because of their mutual proximity and their strengths, 
it is important to note that the correspondence of their 
strengths to the stellar parameters is not straight forward. 

Differences in the carbon spectrum are obtained (De 
Marco et al. 1999b) not only between CMFGEN and 
ISA-wind, but also with the co-moving frame code of 
Koesterke & Hamann (1995). Such discrepancies are not 
due to differences in the atomic data, and it is extremely 
difficult to identify the exact cause, since the effect is 
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Table 2. Summary of model atom used in CMFGEN and ISA-wind radiative transfer calculations, including full levels 
(JVp), super levels (iVs) and total number of transitions (iVrrans), following Hillier & Miller (1998) and de Koter et al. 
(1997; 31 auto-ionizing levels are included for Cm in the atomic data used by ISA- wind) 



Ion iVs JVf A^Ti-ans Details 





CMFGEN 


CMFGEN 


ISA- wind 


CMFGEN 


ISA- wind 


CMFGEN 


ISA- wind 


He 1 


27 


39 


51 


905 


588 


n <14 


n <20 


Hell 


13 


30 


20 


435 


190 


n <30. 


n <20 


ClI 


39 


88 


7 


791 


10 


nl <2s2p3p 2 D 


nl <2s ( S)3d D 
nl <2s( 2 S)5f 1 F° 


Cm 


99 


243 


44 


5513 


174 


nl <2sl0z 1 Z 


Civ 


49 


64 


15 


1446 


105 


n <30 


n <10 


On 


23 


75 




3 




nl <2p 32 P° 




Om 


50 


50 


3 


213 




nl <2p3d 1 D° 


2s 2 2p 2 ( 1 S) 


Oiv 


30 


72 


22 


67 


52 


nl <2s2p4d 2 P° 


2s2p( 3 P°)3d 4 D 

2p( 2 P)3p 1 S 

ls 2 ( 1 S)10d 2 D 


Ov 


31 


91 


24 


835 


57 


nl <2s4f 1 F° 


Ovi 


13 


19 


26 


749 


153 


nl <8d 2 D 


Si iv 


12 


20 




72 




n < 20 




Cam 


16 


39 




104 




nl <3s 2 3p 5 5s 3 P° 




Caiv 


8 


26 




12 




nl <3s 2 3p 4 ( 1 D)3d 2 D 
nl <3s3p 4 ( 4 P)3d 5 D 




Cav 


19 


32 




44 






Cavi 


16 


30 




61 




nl <3s3p 3 ( 3 D)3d 4 G° 
nl <3d 4 ( 1 G°)4p 2 P° 
nl <Z& i { 2 'D )Av 1 ¥ 
nl <3d 2 ( 1 S)4p 2 P° 




Feiv 


21 


280 




3767 






Fev 


19 


182 




1921 






Fe vi 


10 


80 




617 







highly dependent on the parameter space investigated. 
Further investigations are under-way. 

A comparison of the CMFGEN synthetic spectrum with 
the rectified ISO spectrum of 7 Velorum is shown in 
Fig. H Our WR+O model spectrum fits the emission lines 
to a good level of accuracy, demonstrating that line fits 
are consistent also in the IR. Spectral features in the 
H and K windows agree qualitatively with unpublished 
AAT spectroscopy of L.F. Smith (priv. com.). As men- 
tioned in Sect.|[ the combined O+WR model continuum 
is ~10% lower than the flux-calibrated ISO spectrum for 
the adopted light ratio, reddening and observed v magni- 
tude. 

In conclusion it was felt that given the limitations of 
the observed spectra, the parameters of WR11 from the 
optical fits have been constrained rather well. The uncer- 
tainties listed in Table || are determined from the fits (for 
the T*, M, C/O and Voo parameters) or from error propa- 
gation. In the case of the spectroscopic luminosity, L, the 
uncertainty derives directly from the uncertainty on My 
(see Sect. ||) since the uncertainty on the bolometric cor- 
rection BC is negligible. For the radius, i?* , the error prop- 
agates from the luminosity and temperature uncertainties. 
Luminosity, mass-loss and terminal velocity determine the 
uncertainty on the wind efficiency rj. Finally the uncer- 
tainty on the ionizing fluxes, Qo and Qi derive directly 
from L, without taking into account internal model un- 
certainties. 

Given the different nature of the two non-LTE codes, 
it is remarkable how similar the deduced parameters are. 



Table 3. Comparison of stellar parameters obtained for 
WR11 using CMFGEN and ISA- wind (blanketed including 
photon loss) 



Parameters 


CMFGEN 


ISA-wind 


Error 


Filling factor, / 


10% 


100% 




Photon Loss 




yes 




T, [kK] 


57.1 


56.0 


3% 


log(M/v7 [Moyr" 1 ]) 


-4.53 


-4.48 


20% 


C/He[by number] 


0.15 


0.06 


40% 


0/He[by number]' 1 


0.03 


0.01 




R [R©] 


3.2 


3.2 


20% 


log(L,/L ) 


5.00 


4.95 


30% 


BC [mag] 


-4.00 


-3.89 


negl. 


Voo [km s _1 ] 


1550 


1550 


10% 


v/Vf 


21 


28 


60% 


logtQo/s" 1 ) 


48.81 


48.75 


30% 


logtQx/s- 1 ) 


47.76 


47.83 


30% 



"This value follows from the adopted O/C=0.2 



This confirms that despite the approximations adopted by 
ISA-wind, this faster code is suitable for the analysis of 
Wolf-Rayet stars in the WC8 parameter domain (as it was 
also shown to be appropriate for a WN8 star by Crowther 
et al. 1999). 

5.3. UV ionizing fluxes and bolometric luminosity 

In Fig. we present the predicted extreme UV fluxes from 
cmfgen (solid line, top panel) and isa-wind (solid line, 
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Fig. 7. Comparison between the output flux of CMFGEN 
(top) and is A- wind (bottom). Both fluxes are compared to 
the un-blanketed ISA- wind spectrum for the same stellar 
parameters 

bottom panel) compared with the un-blanketed ISA- wind 
model flux (dotted lines) on the energy vs log (AFa) plane. 

As can be seen, the isa-wind flux at A < 400 A is 
harder than the CMFGEN flux, a discrepancy that was al- 
ready pointed out by Crowther et al. (1999). On the other 
hand the Rydberg vs. XF\ plot is designed to emphasise 
the 0-912 A region (1 Rydberg = 13.6 eV) and we see 
by looking at the values of log(Qo/s _1 ) and log(Qi/s _1 ) 
shown in Table that in fact the two fluxes are not par- 
ticularly dissimilar short of 912 A and 504 A. On the other 
hand, the fluxes shortward of ^400 A, which influence the 
strengths of nebular lines such as 2+ , Ar 3+ and Ne 2+ , 
are very different and this discrepancy could affect derived 
nebular properties. 

The spectroscopic luminosity determined by CMFGEN 
and ISA-wind are not dissimilar, although the former is 
slightly higher (Table ||) . Both the luminosities are lower 
than that determined in Paper I from alternative meth- 



ods (log(L/L ) = 5.2±0.1), even considering the formal 
uncertainty determination. 

The luminosity of 63000 L Q obtained by Schaerer et 
al. (1997) can be compared to that obtained here by the 
cmfgen code of 100000 Lq. We can accredit the increase 
to the presence of line blanketing (which includes photon 
loss from the He 11 Lya - see next section), although the 
difference of the two codes and the neglect of carbon and 
oxygen might also play a role. 

6. The effect of line-line interaction 

Line- line interaction between the very strong He II Lya 
and nearby weaker lines has a large effect on the wind ion- 
ization structure. Approximations in the line-blanketing 
technique that affect either the position or the strength 
of even the weakest lines can affect greatly the derived 
stellar parameters. In this section we determine the size 
of this effect and show its consequences for the derived 
parameters of 7 Velorum. 

Schmutz (1997) proposed that capturing of Hen Lya 
line-photons by resonant metal transitions, would reduce 
the ionization balance of the wind. Including line-line in- 
teraction between He 11 Lya and lines from ions such as 
CaV, Fevi, Nivi and Om in his non-LTE model of the 
WN4 star WR6, led to a significantly lower wind ioniza- 
tion structure than that obtained without taking into ac- 
count such interaction. As a result, he had to adopt a 
higher effective temperature to fit the spectrum of this 
star, leading to a larger luminosity and a lower wind per- 
formance number. In this way enough radiation force was 
obtained to drive the outer parts of the wind of WR6. 

Schmutz (1997), however, carried out only an approx- 
imate calculation, which kept the fraction of photons cap- 
tured by metal lines constant throughout the wind, in- 
stead of varying with radius. Additionally, the populations 
of the levels responsible for the metal transitions intercept- 
ing the He 11 photons, were calculated in LTE, instead of 
being calculated consistently with the populations of all 
other levels via the statistical equilibrium equations. Since 
the extent of this effect could lead to the long-sought an- 
swer to a driving mechanism for WR stars, we thought it 
important to investigate it further. 

The line-blanketed cmfgen code (Hillier & Miller 
1999) automatically takes into account line-line interac- 
tion between all the lines included in the adopted atomic 
model so that photon loss should naturally take place in 
our calculation of Sect. between He 11 Lya and the in- 
cluded lines of iron, oxygen, calcium and carbon. To de- 
termine the effect of line-line interaction on the wind ion- 
ization, we removed from the atomic model sets of lines 
in the spectral region neighboring (±2000 km s _1 ) Hen 
Lya. Throughout the rest of this discussion, we will refer 
to the model calculated with this reduced atomic model the 
"no-PL model" (which stands for no photon loss), while 
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the model calculated in Sect. which fits the spectrum of 
7 Velorum, is referred to as the "basic model". 

The elimination of 21 Fevi lines, results in a small in- 
crease in the wind ionization balance with respect to the 
model of Sect.|[ Next we eliminated the O ill Bowen lines 
at 305.72 and 303.65 A ( 3 P- 3 P° and 3 P~ 3 D°) and a fur- 
ther increase in the ionization balance was observed. Next 
we eliminated 1 CaV line, with no change in the calculated 
spectrum and finally we removed 63 C ill lines, which pro- 
duced a further small increase in the ionization balance. In 
Fig. ^| we compare the synthetic spectrum from the basic 
model (solid line) with the spectrum from the no-PL model 
(dashed line), while in Fig. ^ we present the corresponding 
ionization structures. The basic model presents a less ion- 
ized spectrum, with weaker Hen lines and stronger Hei 
lines. The ionization structure of helium is clearly higher 
in the model without the metal lines, although that of car- 
bon and oxygen appear similar, even if the lines of ClV 
follow the trend of the He 11 lines. Eliminating some metal 
lines docs indeed change the wind ionization balance. The 
ionization structure and synthetic spectrum of the basic 
model can be matched of we increase the effective temper- 
ature of the no-PL model by about 10 000 K. 

As a further test of the effect of line-line interaction on 
the overall wind ionization we calculated another model, 
with the same atomic data as in SectJ^, but with the tur- 
bulent velocity of the lines in the He 11 Lya (±300 km s _1 ) 
region set at 5 km s _1 Q instead of 50 km s _1 . In this way 
line-line interaction is reduced, since a smaller turbulent 
velocity leads to narrower lines, smaller overlap and there- 
fore less interaction. We refer to this model as the "low- 
vturb model" 

In Fig. ^ we show the synthetic spectrum for the low- 
vturb model (dotted) compared to the basic model (solid), 
and the no-PL model (dashed), while in Fig. ^ we com- 
pare their respective ionization structures. The ionization 
equilibrium of the low-vturb model is higher than for the 
basic model but similar to the no-PL model. By reducing 
the line widths we have reduced the interaction between 
lines with an overall shift to a higher wind ionization, sim- 
ilar to the shift produced by the elimination of metal lines 
near He II Lya. 

The interception of continuum photons by metal lines 
(continuum blanketing) is an important effect and the 
shift in ionization observed when eliminating metal lines 
(solid and dashed lines in Figs. ^ and ^ could also be 
due to the lack of continuum blanketing by the lines we 
eliminated. The change in ionization between the basic 
model and the low-vturb model (solid and dotted lines in 
Figs. H and ^|), on the other hand, is due to a reduced 
line-line interaction. The similarity of the ionization shifts 
in the two trials (dashed and dotted lines in Figs. || and 
^), indicates that line-line interaction plays an important 
role in the change observed when eliminating metal lines. 

2 This yields a Doppler velocity of 13 km s _1 . 



Weak but numerous metal lines, with wavelengths near 
that of strong emission lines are important in the overall 
ionization balance of the wind and therefore in the derived 
parameters. 

From earlier trials we can estimate that the difference 
between the basic model (solid) and the no-PL or the low- 
vturb models (dashed and dotted, respectively) in Fig. || 
is equivalent to about 10 kK in effective temperature, but 
only about 20% in bolometric luminosity (because the ra- 
dius to which the stellar temperature refers also changes). 
We can therefore conclude that for this parameter space 
photon-loss from the Hen Lya, does lead to a higher lu- 
minosity, although not as high as that determined via the 
evolutionary mass-luminosity relation. 

7. Discussion 

In the following sections we present a critical summary of 
the work carried out in this paper. 

7.1. Stellar parameters and bolometric luminosity 

We have presented a detailed spectroscopic analysis of the 
WR star in the 7 Velorum WR+O binary system. Using 
the wavelength-dependent light ratio and the synthetic O 
star spectrum from Paper I we have generated a corrected 
WR optical and IR spectrum which we have fitted with 
a synthetic spectrum using the CMFGEN non-LTE model 
atmosphere code. Results are summarised in Table U Our 
derived mass-loss is in excellent agreement with the mass- 
loss derived by Stevens et al. (1996) from X-ray ASCA ob- 
servations, once their value is scaled to the hipparcos 
distance. Although the parameters of temperature, mass- 
loss and carbon abundance are well constrained, the oxy- 
gen abundance is assumed. The discrepancy between the 
Hen lines at 4340, 4859 and 6560 A and other Hen lines 
remains to be explained. 

Previous analyses of WR11 were carried out by 
Schaerer et al. (1997) and by De Marco et al. (1999a - us- 
ing the stellar atmosphere code of Koesterke & Hamann, 
1995). Schaerer et al. determined the stellar effective tem- 
perature by fitting the Hei A4471 / Hell A4541 ratio, as- 
suming a (higher) mass-loss and carbon abundance. De 
Marco et al. (1999a) derived larger mass-loss and tem- 
perature, and a lower wind terminal velocity. This can 
be attributed to the fact that in the parameter space of 
7 Velorum the helium spectrum can be reproduced ap- 
proximately with a range of temperature and mass-losses. 
On the other hand, when the lines are analysed in de- 
tail, it is clear that for temperatures larger than about 
57 kK and log(M/M Q yr _1 )>-4.45, predicted line shapes 
deteriorate, becoming wider. This effect was partly com- 
pensated for by De Marco et al. (1999a) by adopting a 
lower wind velocity, although it is clear from their fits 
that most lines are still too wide. Further, their neglect of 
line-blanketing and therefore of photon loss (no line-line 
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Table 4. Summary of the stellar parameters obtained for WR11 using CMFGEN and ISA- wind (blanketed with photon 
loss activated), compared with those derived by Schaerer et al. (1997) and De Marco et al. (1999a) with a non-line 
blanketed model with no photon loss implemented. The results obtained by Dessart et al. (2000) for WR135 are also 
included 



Parameters 


CMFGEN 


ISA-wind 


De Marco et al. 


Schaerer et al. 


WR135: Dessart et al. 


Filling factor, / 


10% 


100% 


100% 


100% 


10% 


Photon Loss 




yes 


no 


no 




T, (kK) 


57.1 


56.0 


75.9 


51.0 


63.0 


log(M/v7) ^©yr" 1 


-4.53 


-4.48 


-4.33 


-4.2 


-4.9 


C/He(by number) 


0.15 


0.06 


0.143 


0.25 


0.13 


0/He(by number) 


0.03 


0.01 






0.03 


R (He) 


3.2 


3.2 


1.84 


3.3 


3.4 


Iog(L,/L Q ) 


5.00 


4.95 


5.01 


4.8 


5.20 


BC (mag) 


-4.00 


-3.89 


-4.1 


-3.5 


-4.0 


Vac (km s -1 ) 


1550 


1550 


1300 


1450 


1400 


v/VT 


21 


28 


29 


144 


5.4 


*For AM V = 1.36. 



interaction between the key lines), might have contributed 
to the different parameters. 

With respect to the analysis of Schaerer et al. (1997), 
the star becomes twice as luminous, 10% hotter and shows 
a mass-loss lower by a factor of two. The carbon abun- 
dance is lower by 40% (although we remind the reader 
that Schaerer et al. assumed both the mass-loss and the 
carbon abundance). This carries implications for studies 
of wind- wind interaction. From this study and the results 
presented in Paper I for the O star, we determine: 



(M v 



oc)WR 



(M Voo)o 



33 



compared to a value of 208 derived by Schaerer et al. 
(1997). Adopting the relationships of Eichler & Usov 
(1993) we determine that the region of stellar wind col- 
lision is approximately 9-30 R© from the O star (0.15 of 
the distance that separates the two stars (63-200 R© - 
Schmutz et al. 1997) vs. 0.06 for the study of Schaerer 
et al. (1997)). Overall the stellar parameters derived are 
not dissimilar from those derived for the single WC8 star 
WR135 by Dessart et al. (2000; see Table |). This indi- 
cates that it is unlikely that substantial levels of emission 
are produced by the collision region or that an anomalous 
ionization is induced by the O star ionizing flux. 

The derived spectroscopic luminosity (log L/Lq = 
5.0 ± 0.1) is increased with respect to that found in pre- 
vious studies, partly due to line-line interaction, and is 
close to the luminosity derived in Paper I via the mass- 
luminosity relationship for WR stars (logL/L© = 5.2 ± 
0.1). An important implication of the higher luminosity 
and of a 10% clumping factor, is that it helps to reduce 
the performance number by a factor of ~20, from 144 to 
7, assuming a 10% clumping filling factor. Lucy & Ab- 
bott (1993) showed that radiation pressure on spectral 
lines should in principle be able to drive winds with per- 



formance numbers ~ 10. So, the strongly reduced 77 sug- 
gests that line driving is the mechanism responsible for 
the dense wind of Wolf-Rayet component in 7 Velorum. 

7.2. CMFGEN VS. ISA- WIND 

Parameters derived with the fast Sobolev approximation 
code ISA- WIND are similar to those derived with CMFGEN, 
apart from the derived carbon mass fraction, with CM- 
FGEN resulting in a carbon abundance which is a factor 
of two larger than ISA- wind. This discrepancy is not to- 
tally understood although we should note that the co- 
moving frame code of Koesterke & Hamann (1995) com- 
pares favourably with CMFGEN for cases in common (e.g. 
for WR135 - Dessart et al. 2000). 

A secondary difference is in the fluxes at A<400 A. 
Crowther et al. (1999) tested the far UV model fluxes of 
CMFGEN and ISA-wind by modelling the H 11 region as- 
sociated with the WN8 star WR124, showing that, for 
that system, the ISA-WIND flux was too hard. Also in 
the current analysis we find that the ISA-WIND flux is 
harder for A<400 A: we suspect that this discrepancy is 
due to the lack of photon branching in the Monte Carlo 
code, which would lead to distributing energetic photons 
to longer wavelengths. Although the ionising properties of 
the two synthetic atmospheres (the values of Qq and Qi) 
are similar, the difference in the extreme UV can carry 
implications when model WR fluxes are used to interpret 
nebular line emission from extra-galactic stellar popula- 
tions (Leitherer et al. 1999). 

We stress that the comparison between the two codes 
is reasonable only if photon loss is included in ISA-WIND 



via the approximation detailed in Sect. 4.3. If photon loss 



were not accounted for, the determined ISA- wind effective 
temperature would be about 8000 K lower, comparable to 
the difference found in SectJ^ when reducing the interac- 
tion between metal lines and He uLya. 
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7.3. The photon loss effect 

The photon loss effect demonstrated by Schmutz (1997) 
as the cause of a lower wind ionization, was tested here 
using the CMFGEN code which treats line-line interaction 
as part of the blanketing. For this parameter space the in- 
teraction between the He II Lya and nearby O in and C ill 
lines causes a shift to lower ionization equilibrium, equiv- 
alent to ~10 000 K or a ~20% increase in the luminosity^. 
By reducing the turbulent velocity from 50 to 5 km s _1 for 
lines residing near 303 A, we have shown that the wind ion- 
ization is enhanced in a similar way, demonstrating that 
line-line interaction, and not continuum blanketing, is at 
the origin of the ionization shift observed when eliminat- 
ing metal lines near the Hen Lya. If line-line interaction 
between strong resonance lines and weaker but more nu- 
merous metal lines leads to a shift in the wind ionization, a 
detailed treatment of even the weakest lines is important. 
We remind the reader that our test remains incomplete 
because, although CMFGEN treats line-line interaction not 
all lines are accounted for. Amongst the ones that are not 
treated there might be some that could play an impor- 
tant role in this mechanism. Additionally, our conclusion 
is only valid within the parameter space appropriate for 
WR11, a WC8 star, since for different temperature and 
mass- loss the strength of the He 11 and other strong lines 
will be different (cf. Crowther et al. (1999) who found that 
photon-loss did not play a role in the atmosphere of their 
WN8 star). 

One should conclude from the presence of the photon 
loss mechanism that a proper treatment of the metal lines 
not only requires one to include very many lines (to real- 
istically model the blanketing), but also the very critical 
ones, such as those near the strong Hen Lya and possi- 
bly at other resonance lines such as Civ at 1550 A. Fur- 
thermore, the presence of photon loss shows that applying 
mean blanketing factors (Pauldrach et al. 1996) - averaged 
over 20-50 A intervals - may miss out some interaction 
with an overall effect on the derived stellar parameters for 
WR stars. Finally, it is clearly pointed out that macro- 
turbulence is connected with photon loss, as the amount 
of turbulence essentially dictates the size and strength of 
the set of lines involved in the line-line interactions. 
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Fig. 4. Comparison between the rectified optical (heros) spectrum of WR11 (solid line) with the synthetic spectrum 
produced by CMFGEN (upper spectrum, dotted line) and ISA-wind (accounting for line blanketing - lower spectrum, 
dotted line). For clarity of display, the regions 4600-4700 A and 5650-5850 A are multiplied by 0.2 and 0.3, respectively, 
and offset 
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Fig. 5. Comparison between WR135 (solid line) and WR11 (dotted line), showing the similarity of the two spectra 
except for the He 11 lines at 4859 and 6560 A 
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Fig. 6. Comparison of the rectified ISO observations of 7 Velorum (solid) with the CMFGEN model (dotted) scaled to 
the optical L A4 7oo(0) / L A4 7oo(WR) = 3.8 
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Fig. 8. Comparison between the model presented in Sect. |5| (solid), the model with metal lines in the Hen Lya region 
removed (dashed) and the model with the reduced vtnrb (dotted) 
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Fig. 9. Comparison between the ionization structure of helium, carbon, oxygen and the level population of He 11 for 
the models presented in Sect. || (solid), the model with metal lines in the Hen Lya region removed (dashed) and the 
model with the reduced v tU rb (dotted) 



